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1. Introduction 

The surface of a sandpile cloaks within it a vast array of complex structures - networks 
of grains whose stability is interconnected, surrounded by pores and necks of void space. 
Bridges - arch-like structures, where mutual stabilisation is a principal ingredient - are 
prime among these, spanning all manner of shapes and sizes through a granular medium. 
They can be stable for arbitrarily long times, since the Brownian motion that would 
dissolve them away in a liquid is absent in sandpiles - grains are simply too large for 
the ambient temperature to have any effect. As a result, they can affect the ensuing 
dynamics of the sandpile; a major mechanism of compaction is the gradual collapse of 
long-lived bridges in weakly vibrated granular media, resulting in the disappearance of 
the voids that were earlier enclosed pQ. Bridges are also responsible for the 'jamming' [2] 
that occurs, for example, as grains flow out of a hopper. 

The difficulty of even identifying, leave alone analysing, structures as complex 
as bridges in a three-dimensional assembly should not be underestimated; such an 
algorithm now exists, and its use has resulted in the identification and classification of 
a panoply of bridge configurations generated via numerical simulations 011]. However, 
even given the presence of such data, it is a far from obvious task to cast it in theoretical 
terms. In the following, we present a theory for the formation of bridges, which, among 
other things, is able to explain at least some of the salient features of our numerical 
data. 

We first define a bridge. Consider a stable packing of hard spheres under gravity, 
in three dimensions. Each particle typically rests on three others, which stabilise it, in 
the sense that downward motion is impeded. A bridge is a configuration of particles 
in which the three-point stability conditions of two or more particles are linked; that is, 
two or more particles are mutually stabilised. They are thus structures which cannot be 
formed by the sequential placement of individual particles; they are, however, frequently 
formed by natural processes such as shaking and pouring, where cooperative effects 
arise naturally. In a typical packing, upto 70 percent of particles are involved in bridge 
configurations. The detailed history of a stabilisation process identifies a unique set of 
bridges, since the stabilising triplet for each sphere is well defined. While it is impossible 
to determine bridge distributions uniquely from an array of coordinates representing 
particle positions, we are able via our algorithm to obtain ones that are the most likely 
to be the result of a given stabilisation process [3J Ej • 

We now distinguish between linear and complex bridges via a comparison of 
Figures [U and [2j Figure [T] illustrates a complex bridge, i.e., a mutually stabilised cluster 
of five particles (shown in green), where the stability is provided by six stable base 
particles (shown in blue). Of course the whole is embedded in a stable network of grains 
within the sandpile. Also shown is the network of contacts for the particles in the bridge: 
we see clearly that three of the particles each have two mutual stabilisations. Figure El 
illustrates a seven particle linear bridge with nine base particles. This is an example of 
a linear bridge. The contact network shows that this bridge has a simpler topology than 
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that in Figure ^ Here, all of the mutually stabilised particles are in sequence, as in a 
string. A linear bridge made of n particles therefore always rests on = n + 2 base 
particles, whereas the number rib of base particles of a complex bridge of size n varies 
from one bridge to another, and always obeys n h < n + 2, because of the presence of 
loops in the contact network of complex bridges. 




Figure 1. A five particle complex bridge, with six base particles (left), and the 
corresponding contact network (right). Thus n = 5 and n\, = 6 < 5 + 2. 




Figure 2. A seven particle linear bridge with nine base particles (left), and the 
corresponding contact network (right). Thus n = 7 and n\, = 9 = 7 + 2. 

An important point to note is that bridges can only be formed sustainably in 
the presence of friction; the mutual stabilisations needed would be unstable otherwise! 
Although we use Monte Carlo simulations (described below) which do not contain 
friction explicitly, the configurations we generate correspond to those in nature that 



Statistics of bridge geometries 



4 



include frictional effects. In particular, they generate coordination numbers in a range 
that is consistent with the presence of friction P3 |SJ E] . Additionally, we have analysed 
the configurations of molecular dynamics simulations in the limit of high friction§; these 
generate distributions of bridges very similar to our own. Our simulations suggest strong 
analogies between bridge structure and force distributions in granular media, which have 
been measured both in experiments [7j and in simulations [S] . This indicates that bridges 
might be complementary (and experimentally easier) probes of inhomogeneities in these 
systems. 

2. Simulation details 

We have examined bridge structures in hard assemblies that are generated by an 
established, non-sequential, restructuring algorithm [TJ [HI E] , whose main modelling 
ingredients involve stochastic grain displacements and collective relaxation from them. 
This algorithm restructures a stable hard sphere deposit in three distinct stages. 

(1) The assembly is dilated in a vertical direction (with free volume being introduced 
homogeneously throughout the system), and each particle is given a random 
horizontal displacement; this models the dilation phase of a vibrated granular 
medium. 

(2) The packing is compressed in a uniaxial external field representing gravity, using a 
low-temperature Monte Carlo process. 

(3) The spheres are stabilised using a steepest descent 'drop and roll' dynamics to find 
a local minimum of the potential energy. 

Steps (2) and (3) model the quench phase of the vibration, where particles relax to 
locally stable positions in the presence of gravity. Crucially, during the third phase, 
the spheres are able to roll in contact with others; mutual stabilisations are thus 
allowed to arise, mimicking collective effects. The final configuration has a well-defined 
network of contacts and each sphere is supported, in its locally stable position, through 
point contacts, by a set of three other spheres uniquely defined. In practice, the final 
configuration may include a few non-stabilised particles. 

The simulation method recalled above builds a sequence of static packings. Each 
new packing is built from its predecessor by a random process and the sequence achieves 
a steady state. In the steady state thus obtained, structural descriptors such as the mean 
packing fraction and the mean coordination number fluctuate about well-defined mean 
values. The steady-state mean volume fraction $ is found to lie typically in the range 
$ ~ 0.55 - 0.61. This volume fraction depends on all the parameters of the simulation, 
and can therefore be adjusted to any desired value in the above range. The mean 
coordination number is always Z rs 4.6 ±0.1. Since for frictional packings the minimal 

§ We thank Leo Silbert for letting us use configurations generated from his molecular dynamics 
simulations, where we have found and analysed the distribution of bridges. 
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coordination number is Z = d+1 [5], this confirms that our 3D configurations correspond 
to those generated in the presence of friction. In fact, the mean coordination number 
of molecular dynamics configurations of frictional sphere packings is slightly above 4.5, 
in the limit of a large friction coefficient We recall that, in the absence of friction, 
the coordination number is that of an isostatic system, Z = 2d P|, hence Z = 6 in the 
present three-dimensional situation. 

To be more specific, simulations were performed in a rectangular cell with lateral 
periodic boundaries, and a hard disordered base. The simulations, performed serially 
on a desktop workstation, have a very time-consuming Monte Carlo compression phase 
with each simulation taking several days' worth of CPU time. This is necessary in order 
to have simulation results which are reproducible, without appreciable dependence on 
Monte Carlo parameters or system size. Each of our configurations includes about 
iVtot = 2200 particles. We have examined approximately 100 configurations from the 
steady states of the reorganisation process for the following two values of the packing 
fraction: <3> = 0.56 and <£> = 0.58. Segregation is avoided by choosing monodisperse 
particles; a rough base prevents ordering. A large number of restructuring cycles 
is needed to reach the steady state for a given shaking amplitude; about 100 stable 
configurations (picked every 100 cycles in order to avoid correlation effects) are saved 
for future analysis. From these configurations, and following specific prescriptions, our 
algorithm identifies bridges as clusters of mutually stabilised particles j31 Ej ■ 

Figure El illustrates two characteristic descriptors of bridges used in this work. Along 
the lines of [3], the main axis of a bridge is defined using a triangulation of its base 
particles. Triangles are constructed by choosing all possible connected triplets of base 
particles: the vector sum of their normals is defined to be the direction of the main axis 
of the bridge. The orientation angle B is defined as the angle between the main axis 
and the z-axis. The base extension b is defined as the radius of gyration of the base 
particles about the z-axis (which is thus distinct from the radius of gyration about the 
main axis of the bridge). 

3. Geometrical characteristics. Size and diameter distribution 

In the following, we present statistics for both linear and complex bridges; however it is 
the former that we analyse in greater detail, both because they are conceptually simpler 
and because they are in fact more numerous in our simulations. The latter clearly 
show that as a linear bridge gets longer, it is increasingly likely to develop branches, 
whereupon it becomes complex [3I31- 

The formation of either type of bridge in shaken granular media is, of course, a 
dynamical process; we, however, adopt an ergodic viewpoint [Ej here, inspired by the 
ensemble statistics of bridges generated by our simulations. As in polymer theory [llj, 
we differentiate between linear and complex geometries; in our phenomenology, we 
regard a linear bridge as a random chain which grows in arc length s by the sequential 
addition of links to its end, while complex bridges, being branched structures, cannot 
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Figure 3. Definition of the angle and the base extension b of a bridge. The 
main axis makes an angle with the z-axis; the base extension b is the radius 
of gyration of the base particles about the z-axis. 



grow sequentially. This replacement of what is in reality a collective phenomenon by a 
random walk is somewhat analogous to the 'tube model' of linear polymers ^T]: both 
are simple but efficient effective pictures of very complex problems. 

In this spirit, we first address the question of the length distribution of linear 
bridges. We define the length distribution f n as the probability that a linear bridge 
consists of exactly n spheres. We make the simplest and the most natural assumption 
that a bridge of size n remains linear with some probability p < 1 if an (n + l) th sphere 
is added to it, in the above sense. This simple effective picture therefore leads to the 
exponential distribution 

f n = (l-p)p n . (3.1) 

This exponential law is a common feature of all models for random 'strings' made of n 
units, irrespective of the way they are formed. An illuminating example is provided 
by percolation clusters in one dimension, which are the prototype of random linear 
objects [12]. 

The above argument leading to an exponential distribution can alternatively be 
reformulated as follows, by means of a continuum approach. This will give us the 
opportunity to introduce a formalism that will be extensively used later. A linear bridge 
is now viewed as a continuous random curve or 'string', parametrised by the arc length s 
from one of its endpoints. Within the above picture of a sequential process, we assume 
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that a linear bridge disappears at a constant rate a per unit length, either by changing 
from linear to complex or by collapsing. The probability S(s) that a given bridge 
survives at least up to length s thus obeys the rate equation S = —aS, and therefore 
falls off exponentially, according to S(s) = exp(— as). The probability distribution of 
the length s of linear bridges therefore reads f(s) = —S(s) = a exp(— as). This is the 
continuum analogue of (j3.1|) . 

Figure 0] shows a logarithmic plot of numerical data for the length distribution f n of 
linear bridges. The data exhibit an exponential fall off of the form (|3.1j) . with p »s 0.37, 
i.e., a ~ 0.99||. This exponential decay of the distribution of linear bridges is very 
accurately observed until n « 12. 




n 

Figure 4. Logarithmic plot of the length distribution /„ of linear bridges, 
against n, for <I> = 0.58. Full line: least-square fit yielding p ~ 0.37, i.e., 
a » 0.99. 



Complex bridges begin to be dominant around n « 8. For a complex bridge, 
the number n is now referred to as the 'size' of the bridge. As complex bridges 
are branched objects, it is natural to expect that their size distribution resembles 
that of generic randomly branched objects, such as e.g. lattice animals or critical 
percolation clusters ^2] • These prototypical examples have a power-law fall-off in their 
size distribution: 

fn ~ n T . (3.2) 

The numerical data shown in Figure clearly show that the size distribution of bridges 
obeys a power- law behaviour of the form ()3.2j) . and thus behave as generic randomly 

|| Despite its proximity to unity, this value of a has no special significance, as can be seen from the 
corresponding value of p. 
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branched objects. The measured value of the exponent r seems to coincide with the 
limiting value r = 2, at and below which the mean size (n) = n fn is divergent. 
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Figure 5. Log-log plot of the size distribution f n of all bridges, against n, for 
two volume fractions. This distribution is dominated by complex bridges, since 
these predominate after n ~ 8. Full line: least-square fit yielding r ~ 2.00. 



We now turn to the diameter of linear and complex bridges, which is one of the 
most important of their geometrical characteristics. The typical diameter R n of a bridge 
of size n is such that R 2 n is the mean squared end-to-end distance over all the bridges of 
size n. We examine this quantity separately for linear bridges, as well as the ensemble 
of all bridges. In both cases we expect and find a power-law behaviour 

R n ~ n u . (3.3) 

For linear bridges, D\ in = 1/ v u n can be interpreted as their fractal dimension. The 
observed value (see Figure EJ) v\ m « 0.66, i.e., Am ~ 1.51, lies between those for a 
self-avoiding walk in two (u 2 = 3/4) and three (z/ 3 m 0.59) dimensions [T^]. This seems 
entirely reasonable, since linear bridges are likely to behave as typical self-avoiding 
curves. Furthermore, in three dimensions, their effective dimensionality is expected to 
lie between 2 and 3: they may start off confined to a plane, and then collapse onto each 
other, due to the effects of vibration. 

For complex bridges, we again find a power law of the form ()3.3|) . The geometrical 
exponent z/ a u « 0.74 (see Figure EJ) is now naturally compared with the value v ~ 0.875 
of 3D percolation |12j . (Although the figure shows statistics for 'all' bridges, these 
mostly reflect the behaviour of complex bridges, which dominate beyond about n ~ 8.) 

Lastly, we plot in Figure |H1 the probability distribution function of the base 
extension b, for linear bridges; this is clearly a measure of the spanning, and hence 




ln(n-l) 

Figure 6. Log-log plot of the typical diameter R n of linear bridges, against n— 1, 
the number of mutually stabilising bonds, for <3? = 0.58. Full line: least-square 
fit yielding v\\ n sa 0.66, i.e., D\\ a m 1.51. 




12 3 4 5 

In n 

Figure 7. Log- log plot of the typical diameter R n of all bridges, against size 
n, for $ = 0.58. These statistics mainly derive from complex bridges, which 
dominate beyond about n ~ 8. Full line: least-square fit to points with n > 5, 
yielding w 0.74. 

the jamming, potential of a bridge (see Figure EJ. We plot also in Figure |H1 distribution 
functions that are conditional on the bridge size n. This was done to compare our results 



Statistics of bridge geometries 



10 



with those of fn], bearing in mind of course that our results refer to sphere packings in 
three dimensions whereas theirs refer to disk packings in two. In our case, the conditional 
distributions are sharply peaked, and only extend over a finite range of b; in theirs the 
peaks are smaller and the distributions broader. Our results indicate that at least in 
three dimensions, and upto the lengths of linear bridges probed, bridges of a given 
length n have a fairly characteristic horizontal extension; we can predict fairly reliably 
the dimension of the orifice that they might be expected to jam. The results of ^3] 
indicate, on the contrary, that there is a wide range of orifices which would be jammed 
by disk packings of a given length, in two dimensions. Also, for our three-dimensional 
bridges, the cumulative distribution has a long tail at large extensions, reflecting chiefly 
the existence of larger bridges than in the data of [Tlj. (We note that this long tail 
is characteristic of three-dimensional experiments on force chains which we will 
discuss at greater length below.) All of this seems entirely reasonable given the more 
complex topology of three-dimensional space, as well as the greater diversity of possible 
sphere packings within it. 

In the lower part of Figure |H1 we have plotted the logarithm of the normalised 
probability distribution of base extensions, against the normalised variable b/{b), 
where (6) is the mean extension of bridge bases. This figure emphasises the exponential 
tail of the distribution function, and also shows that bridges with small base extensions 
are unfavoured. It is interesting to note that the form of the normalised distribution 
resembles that found in force distributions of emulsions JH]; and that there are also 
similarities with normal force distributions in MD simulations of particle packings 
(see e.g. jHJ [T3| ) . In our simulations, the sharp drop at the origin with a very 
tangible peak strongly resembles that obtained in the latter case, in the limit of strong 
deformations Realising that the measured forces jHl El propagate through chains 
of particles, we use this similarity to suggest that our bridges are really just long-lived 
force chains, which have survived despite strong deformations. We suggest also that 
with the current availability of 3D visualisation techniques such as NMR ^7j, bridge 
configurations might be an easily measurable, and effective, tool to probe inhomogeneous 
force networks in shaken sand. 

4. Orientational distribution of linear bridges. Theory and simulation 

Linear bridges predominate at small n, which is after all the regime of interest for 
finite-sized simulations and experiments. The theory we present below focuses on such 
bridges; the distribution of their orientation is predicted, and then compared with our 
simulation results. 

As before, a linear bridge is modelled as a continuous curve, parametrised by the 
arc length s from one of its endpoints. We propose the following effective description 
of the local and global orientation of this curve. For the sake of simplicity, we choose 
to focus on the most important degree of freedom, which is the tilt with respect to 
the horizontal; the azimuthal degree of freedom is therefore neglected. Accordingly, we 
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Figure 8. Distribution of base extensions of bridges, for $ = 0.58. In the 
top panel, the distributions conditional on the bridge size n, for n = 2, 4, 6 
are marked in blue, while the cumulative distribution is drawn in black. The 
bottom panel shows the logarithm of the normalised probability distribution, 
as a function of the normalised variable b/(b). 



define the local or 'link' angle 8(s) between the direction of the tangent to the bridge at 
point s and the horizontal, and the mean angle made by the bridge from its origin up 
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to point s, 

9(s) = - [ S 0(u)du. (4.1) 
s Jo 

The local angle 9(s) so denned may be either positive and negative; it can even change 
sign along the random curve which represents a linear bridge. Of course, the orientation 
angle measured in our numerical simulations pertains to the three-dimensional realm 
of real sphere packings. It is there positive by construction, being defined at the end of 
Section 2 as the angle between the main bridge axis and the z-axis. 

Our simulations show that the angle 0(s) typically becomes smaller and smaller as 
the length s of the bridge increases. Small linear bridges are almost never flat |31 E] ; as 
they get longer, assuming that they still stay linear, they get 'weighed down', arching 
over as at the mouth of a hopper |18j . Thus, in addition to our earlier claim that long 
linear bridges are rare, we claim further here that (if and) when they exist, they typically 
have flat bases, becoming 'domes'. 

We use these insights to write down equations to investigate the angular distribution 
of linear bridges. These couple the evolution of the local angle 8(s) with local density 
fluctuations 0(s) at point s: 

= -a0-ty 2 + Ai77i(s), (4.2) 
0= -c0 + A 2 r ]2 {s). (4.3) 

The effects of vibration on each of 9 and are represented by two independent white 
noises T)x(s), rj 2 (s), such that 

(Vi(s)Vj(s')) = 2 6^6(3-3'), (4.4) 
whereas the parameters a, A 2 are assumed to be constant. 

The phenomenology behind the above equations is the following: the evolution 
of 9(s) is caused by the addition to the bridge of single particles to the bridge, along 
the lines of our effective sequential picture. The motion of particles within their cages 
gives rise to the fluctuations of local density at a point s; thus <fr can be regarded as an 
effective collective coordinate, with 9 seen as an independent-particle coordinate [TH] . 

The individual dynamics are simple: the first terms on the right-hand side 
of 1)4. 2j) . (|4.Hjl say that neither 9 nor is allowed to be arbitrarily large. Their coupling 
via the second term in (J4.2|) arises as follows: large link angles 9 will be increasingly 
unstable in the presence of density fluctuations 2 of large magnitude, which would to 
a first approximation 'weigh the bridge down', i.e., decrease the angle 9 locally. Density 
fluctuations (otherwise known as 'dilatancy') also play a key role in the dynamics 
of the angle of repose of a sandpile, which are described in a concurrent paper [20J. 
There, sandpile collapse is described in terms of an activated process, where an effective 
temperature competes against configurational barriers generated by dilatancy. 

Reasoning as above, we therefore anticipate that for low-intensity vibrations and 
stable bridges, both density fluctuations (f>(s) and link angles 0(s) will be small. 
Accordingly, we linearise ()4.2|) . obtaining thus an Ornstein-Uhlenbeck equation fH\ 12*2*] 

= -a9 + Ait;i(s). (4.5) 
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Let us make the additional assumption that the initial angle 6^, i- e -> that observed 
for very small bridges, is itself Gaussian with variance o~ 2 = (9q)- The angle 8(s) is then 
a Gaussian process with zero mean for any value of the length s. Its correlation function 
can be easily evaluated to be 

(0(s)9(s f )) = a e 2 q e- a| -'' + (a 2 - <)e- a ( s+s '>. (4.6) 

The characteristic length for the decay of orientation correlations is £ = 1/a, and the 
variance of the link angle reads 

<^>(*) = < + K 2 -Oe-^ (4-7) 

where 

A 2 

a 2 cq = -± (4.8) 

is the 'equilibrium' value of the variance of the link angle. Thus as the chain gets longer, 
the variance of the link angle relaxes from a 2 (that for the initial link) to <r 2 q , in the 
limit of an infinite bridge. 

We predict from the above that the mean angle 0(s) will have a Gaussian 
distribution, for any fixed length s. By inserting ()4.6j) into (|4.1j) . we derive its variance: 

< e2 >M = 2 < M "„v e "" + (-o 2 - O 11 ^ 1 - (4-9) 

The asymptotic result 

9rr 2 OA 2 

<0 2 >( s )^f^^i (4.10) 
as a z s 

confirms our earlier statement (see below ()4.1Jl ) that a typical long bridge has a base 
that is almost flat. It can be viewed as consisting of a large number as = s/£ ^> 1 of 
independent 'blobs'^f, each of length £. 

The result (14.10}) has another interpretation. As O(s) is small with high probability 
for a very long bridge, its extension in the vertical direction reads approximately 

Z = z(i)-z(0)«s9(s), (4.11) 

so that {Z 2 ) m s 2 (6 2 )(s) ~ 2(A!/a) 2 s. Going back to the discrete formalism, we have 
therefore 

Z n ~ n 1 / 2 . (4.12) 

The vertical extension of a linear bridge is thus found to grow with the usual random- 
walk exponent 1/2, whereas its horizontal extension exhibits the non-trivial exponent 
fii n ~ 0.66 of Figure |H1 Thus, long linear bridges are domelike; they are vertically 
diffusive but horizontally superdiffusive. We recall that the two-dimensional arches 
found in were diffusive in a vertical direction; our present results show that in three 
dimensions, the vertical diffusivities of bridge structures remain, and are enriched by a 

This is qualitatively reminiscent of a similar problem in polymers, where the de Gennes picture of 
'blob' dynamics in dilute solutions gives way to rigid rod- like behaviour - see e.g. 
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horizontal superdiffusivity. Evidently, jamming in a three-dimensional hopper would be 
caused by the planar projection of such a dome. 

We now compare these predictions with data from our simulations. The numerical 
data shown in Figure |H1 confirm that the mean angle has, to a good approximation, a 
Gaussian distribution in this typical situation of bridges of size n = 4. We recall that the 
angle measured in simulations pertains to the three-dimensional world. Hence only 
the positive half of the Gaussian is used for comparison with results of simulations, and 
the sinB three-dimensional Jacobian is divided out. It has also been checked to high 
accuracy (less than 2 percent for a bin size Sip = 5 degrees) that the azimuthal angle ip 
of the main bridge axis has the expected flat distribution. Figure ITU1 shows the measured 
size dependence of the variance (0 2 )(s), for both volume fractions. The numerical data 
are found to agree well with a common fit to the first (stationary) term of (|4.9J) . with 
a common value of the parameters a^ q = 0.093 and a = 0.55. The 'transient' effects 
of the second term of ()4.9|) are invisible with the present accuracy (see below). This 
observed agreement therefore provides a first confirmation of the validity of our theory, 
quite remarkable given that it is entirely independent of the simulations. We conclude 
that in spite of its simplicity, our theory captures the main features of linear bridges. 
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Figure 9. Plot of the normalised distribution of the mean angle (in radians) 
of linear bridges of size n = 4, for both volume fractions. The sin© Jacobian 
has been duly divided out, explaining thus the larger statistical errors at small 
angles. Full lines: common fit to (half) a Gaussian law. 

We end this section with the following remarks. First, more subtle effects, including 
the effects of transients via the second term of ()4.9j) . and the dependence of the 
parameters a^ q and a on the packing fraction $, are deserving of further investigation. 
As a preface to our second remark, we recall that in our phenomenological picture 
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Figure 10. Plot of the variance of the mean angle of a linear bridge, against 
size n, for both volume fractions. Full line: common fit to the first (stationary) 
term of (|4.9jl . yielding a^ q = 0.093 and a = 0.55. The 'transient' effects of the 
second term of (|4.9|) are invisible with the present accuracy. 



linear bridge formation proceeds sequentially, while complex bridges with their branched 
structures form collectively, in space. + We might expect that with increasing density $, 
branched structures would become more and more common; linear bridge formation, 
with its 'sequential' progressive attachment of independent blobs would then become 
more and more rare. Our theory should therefore cease to hold at a limit packing 
fraction $ 1iT T1 , which is reminiscent of the single-particle relaxation threshold density [23J, 
at which the dynamics of granular compaction crosses over from being single-particle 
(sequential) to collective. Finally, while we have allowed for the local variability of 
density fluctuations (ft via (j4.Hj) . we have so far not explicitly coupled these to bridge 
orientations. All of these issues form part of ongoing work, since much more numerical 
data is needed before a detailed comparison with theory is possible. 

5. Discussion 

In the above, we have classified bridge structures in granular media as either linear 
or complex, analysed their structures as obtained in our simulations, and obtained 
values for some of their more important geometrical exponents. We have also presented 
a theory for their orientational correlations, which is in reasonably good agreement 
with our simulations, and more importantly, makes predictions that can be observed 

+ We reiterate once again that this effective picture replaces the complexity of the full dynamical 
picture, where all bridges form cooperatively. 
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experimentally. Do linear bridges really predominate at short lengths? Do they really 
cross over to being domelike at large lengths, if they survive? The answers to these and 
other questions, if obtained experimentally, will not just satisfy what is arguably mere 
scientific curiosity. Our investigations above suggest that long-lived bridges are natural 
indicators of sustained inhomogeneities in granular systems. Consequently, experimental 
and theoretical explorations of geometry and dynamics in stable bridge networks will 
have immediate implications for issues such as hysteresis and cooperative stability, which 
are some of the most visible manifestations of complexity in granular media. 
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